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Abstract 

In this paper, a unified approach to sequence-based control and estimation of linear networked 
systems with multiple sensors is proposed. Time delays and data losses in the controller-actuator- 
channel are compensated by sending sequences of control inputs. The sequence-based design 
paradigm is further extended to the sensor-controller-channels without increasing the load of the 
network. In this context, we present a recursive solution based on the Hypothesizing Distributed 
Kalman Filter (HKF) that is included in the overall sequence-based controller design. 



1. Introduction 

With advances in the development and distribution of modern network technologies such as 
Ethernet or WLAN (IEEE 802.11), general purpose networks for cordless communication have 
become an easy and cheap alternative to complex field buses. Although new flexibility is acquired, 
the communication is typically less reliable than in wired real-time architectures and, hence, users 
are faced with additional challenges regarding network-specific disturbances such as transmission 
delays and data losses. In control theory, these kinds of network-related problems are investigated 
in the area of Networked Control Systems (NCS). 

In the NCS community, several approaches have been proposed that consider transmission delays 
and stochastic data losses in the controller design (see [l[ for an overview). In this contribution, we 
focus on a design philosophy called sequence-based control, which stems from the idea that modern 
digital communication networks usually transmit data in form of atomic packets, which enforce 
that either all data of a packet is received or none. Therefore, the idea is to not only send the usual 
data over the network, but also extra information that is used to mitigate the network-induced 
effects. However, the transmission of extra information increases the network load and in general 
causes higher transmission delays and/or loss rates. Nevertheless, it is a common assumption that 
transmitting additional information within a data packet leads to an increased system performance 
as long as the additional information is not too extensive. 
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For the connection between controller and actuator (CA-channel), additional information is 
usually attached in form of predicted control inputs that are applicable at future time steps. The 
actuator stores these control inputs in a buffer so that it can fall back upon a predicted control 
input in case a data packet is lost or delayed. 

In literature, three major lines of sequence-based controller designs can be distinguished. The 
first line of methods is based on a nominal controller that is designed for the system without 
consideration of networked-induced effects and that is extended afterwards to generate sequences 0, 
0|. The second line of methods stems from Model Predictive Control (MPC) theory, where the 
control sequences are generated as byproduct of solving an open-loop receding-horizon optimization 
problem [4]. In the third line of methods, the sequence-based NCS is formulated as a stochastic 
optimal control problem that is solved offline for some system classes 0, [f| . 

In this contribution, we extend a result of the third line that has been derived in 0j for linear 
systems and TCP-like network connections in two ways (see Sec. [2] for comments on TCP-like 
networks). First, we extend the system to multiple independent sensors that individually send 
partial state measurements over the network. We show that in this scenario the separation principle 
holds and discuss the structure of the optimal control law that depends on the minimum-mean- 
squared-error estimate of the state. Second, we apply the sequence-based design philosophy of the 
CA-channel to the data connection between sensors and controller (SE-channel). 

An intuitive approach for the latter extension is that each sensor not only sends the most recent 
measurement to the controller, but also a (finite) set of measurements obtained in previous time 
steps 0], This way, the estimation performance is improved since in case a sensor transmission 
has been dropped or delayed by the network, the missing information is also part of following 
transmissions. But, as pointed out above, the sent measurement sequences should not be too large. 
In particular, in the presence of multiple sensors, the network load increases considerably with the 
number of measurements transmitted per packet. 

Hence, in this paper we investigate an approach how information contained in (possibly infinitely 
long) measurement sequences can be comprised in recursive variables. For that reason, we assume 
the sensors to be capable of performing minor processing tasks. While the Kalman Filter is the 
optimal estimator when all measurements can be processed at one node the estimation quality 
is no longer optimal when local Kalman Filters are employed on multiple sensors 0]. 

As we are interested in an estimation principle that is equivalent to processing the respective 
measurement sequences, we utilize recent results in estimation theory, where it has been shown 
in form of the so called Distributed Kalman Filter (DKF) 0, 0] that measurement data can be 
processed by a group of independent local sensors in a way that the fused result still yields a 
globally optimal result. 

Unfortunately, in the presence of lossy communication networks, the DKF does not work as 
it is necessary to have a complete set of all sensor data at the controller side to generate a state 
estimate. A generalization that manages to provide estimates even under uncertain conditions is 
the Hypothesizing Distributed Kalman Filter (HKF) [lol. [ill]. The second part of this contribution 
focuses on the extension and integration of the recursive HKF algorithm into the sequence-based 
controller design. 

The following of the paper is structured as follows. After the problem formulation in Sec. 
we generalize the sequence-based control result from 0] in Sec. E] by showing that the separation 
principle holds for further classes of available measurement information. The main part of this 
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Figure 1: The considered scenario with one actuator and multiple sensors that are connected to the controller via a 
network. 

contribution is presented in Sec. HI where we derive a generalization of the HKF in order to extend 
the sequence-based methodology to the measurement channel. In particular, we propose novel 
recursive sum formulas for the HKF, introduce a method to initialize estimates from measurements, 
and derive an extension of the HKF that is capable of handling the subsequent inclusion of control 
inputs at the controller side. After discussing the application of the proposed algorithm, we give a 
short summary and an outlook in Sec. 



2. Problem Formulation 

In this paper, we consider the system setup depicted in Fig. [T] consisting of a plant that is 
controlled and observed over a digital network. The plant is supposed to be linear and time- 
invariant and its output is not directly observable but measured by M G N>o independent sensors^ 
according to 



x k+l = Ax k + Bu k + w k , 
i k = C { x k + wj. , 



(1) 
(2) 



where x k G 



U.k € 



and y\ G M 9i denote the plant's state, the control input applied by 



the actuator and the measured output of the i-th sensor (with i G N>o and i < M). The matrix 
A G R nxn is assumed to be regular. The matrices A, B G M. nxm , and C* G R* xn are known by 
the controller and the sensors. The terms w k G M n and v k G M Qi represent mutually independent, 
zero-mean, Gaussian white noise processes with covariance matrices S G R nx ™ and 0* G M qiXqi . 
The initial state of the plant is assumed to be Gaussian distributed with mean ~x and covariance 
matrix Pq. 



The setup can easily be extended to the case of multiple actuators and a time-variant plant. In order to give a 
concise description of the unified sequence-based approach we limit ourselves to this simpler setup. 
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The digital network connections are subject to time-varying transmission delays and stochastic 
packet losses. The stochastic characteristics of these effects are assumed to be known. Furthermore, 
a TCP-like protocol is used for transmissions in the CA-channel. This means that successfully 
transmitted data packets to the actuator are acknowledged at the controller within the same time 
For transmissions between sensors and controller we relax this assumption and consider the 
less restrictive case that no acknowledgments are available. Finally, we assume the network nodes 
to be synchronized, data packets to be tagged with a time-stamp, and controller, actuator, and 
sensors to be time-triggered. 

To compensate for time delays and data losses between controller and actuator, at every time 
step the controller sends control sequences to the actuator that not only contain a control input in- 
tended for application in the current time step M/dfe; but also Na S No control inputs for consecutive 
future time instants. Such a control packet is referred to as \J_\~ and is of the form 



Hi 



T T T 
^k\k ^k+l\k ■■■ l±k+N A \k 



(3) 



where, for example, u.k+\\k denotes a control input calculated at time step k that is intended to be 
applied at time step k + 1. The actuator is equipped with a buffer to store the control sequence 
with the most recent information (among all received sequences). By applying the control input 
of this buffered sequence that matches the current time step, time delays and losses of subsequent 
data packets can be compensated until a packet with more recent information arrives. In case the 
actuator does not receive a new control sequence until the buffer runs out of applicable control 
inputs, the actuator applies a known default control input u d . 

As a counterpart to the CA-channel, network effects between sensors and controller are also 
compensated by a sequence-based design philosophy. More precisely, the sensors not only use the 
measurement of the current time step k to generate a data packet but also include measurements 
of the last Ng £ No time steps. The information used to generate the output packet of the i-th 
sensor ^° ut,J can formally be described by 

z° ut ' 1 = g(yl- Ns :k) , (4) 

where g(-) denotes an arbitrary sensor algorithm and the notation y a -b denotes the set of measure- 
ments {yk\a < k <b}. The union of all M sets Z^ ut ' 1 is denoted by 

z out = z ont,l u z out,2 (J _ (J z ont,M _ (5) 

Due to network disturbances, the controller may receive no, one, or even more than one packet per 
sensor at each time step. The set of packets received by the controller of the i-th sensor at time 
step k is defined by the set Z k ' and the union of these sets at time step k over all M sensors by 

^k ■ 



2 A TCP-like network does not reflect real Ethernet-TCP/IP networks, where acknowledgments may be subject 
to a considerable time delay. However, it is possible to realize a TCP-like network connection by, e.g., message 
prioritization. Furthermore, TCP-like networks are also of theoretical interest since they give insights into the far 
more complex cases of real TCP connections or connections where no acknowledgments are provided by the network 
(as, e.g., for Ethernet UDP/IP networks). 
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The problem investigated in this paper is to find an admissible control law that minimizes the 
quadratic cost function 



f K-l 

= Ei XkQx k + XkQx-k + ul R-Mfc 

I k=0 



LLo:K-l,X ,Po \ , (6) 



where K € N>o is the final time step and Q and R are symmetric weighting matrices that are 
positive semi-definite and positive definite, respectively. A control law is called admissible if it 
requires only information that is available to the controller. In our case, the information available 
to the controller is described by the information set 

Zfc = jl , Po, Z'hkiU.Q-.k-l, #0:fc-l} , (7) 

where Ok represents the information provided by the acknowledgments of the TCP-like connection 
between controller and actuator. 



3. Optimal Control Law 

To solve the sequence-based optimal control problem formulated in Sec. we use a result 
derived in 0], where a simpler problem setup was considered containing only one sensor that 
transmits one raw measurement per time step. It turns out that the results can easily be extended 
to the considered setup with multiple sensors that send (possibly pre-processed) measurement 
information. We summarize the extended result in the following theorem that is formulated in 
terms of the augmented state 



(8) 



that contains the state of the plant and all control inputs of already sent control sequences that 
still could be applied to the plant. The last entry is the default control input. 

Theorem 3.1 Consider the problem of finding an optimal control law to generate control sequences 
of length N a & N>o minimizing the cost function ^ subject to the available information 0, system 
dynamics (pTJ) ^ actuator logic described in Sec.\^ and M sensors that process data according to ^ 
and dl]). Then, 

1. similar as in standard LQG control, the separation principle holds, i.e., the optimal control law 
can be separated into a) an optimal state estimator that calculates the conditional expectation 
E \ $,Alk \ and b) an optimal state feedback controller that utilizes the feedback matrix L^, 



y— k\k-l —k+l\k-l 

rp rp 

i%|fe-2 — fc+l|fc-2 



rp rp 
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2. the optimal control law is linear in the conditional expectation of the augmented state according 



3. and the state feedback matrix explicitly depends on the delay probability distribution of 
the controller-actuator-network. It can be calculated as given in Theorem 1 of [B] since it is 
identical to the feedback matrix for the same system with only one sensor sending one raw 
measurement per time step. 

Proof. When we replace the measurement equation and information set If. of 0] with the corre- 
sponding equations and (J7]) of our problem, Lemma 1 of 0] still holds. Therefore, the conditional 
expectation E j^jZfc j is stochastically independent of the control sequence U_ k . This implies that 
the separation principle holds and thus, proves part one of the theorem. Part two and three can 
be proved analogously to the proof of Theorem 1 in with replaced by © . □ 



Remark 3.1 It is interesting to note that Theorem \'J.l\ also holds for all information structures that 
can be characterized by cumulative subsets X k of the information setXk, i.e., by subsets that satisfy 
Ife-l Q X-k an d Ik Q%k- This holds regardless of the choice of Ng, i.e., how much information of 
previous time steps is transmitted per packet. Therefore, Theorem \3.1\ also holds for the case of 
infinite measurement sequences. 

According to Theorem I3.1} the optimal control law consists of the combination of an optimal 
state estimator and an optimal state feedback controller. Since the optimal state feedback controller 
is the same as in [5], we focus on finding the optimal sequence-based state estimator to calculate 
the conditional expectation E|£j2fcj in the following. 

4. Optimal Estimator 

As discussed in the introduction, one way to implement a sequence-based information approach 
for the SE-channel is to include the last N$ measurements into the sensor output packet at every 
time step. This corresponds to setting g = id in and results in 



to 



U-k = Lfc • E 




(9) 





(10) 



and thus, the estimation accuracy of the conditional mean 



E 



{i k \ik}=v{i k \uN S )} 
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is best when the complete sequence of measurements is transmitted to the controller, i.e., Ng = k. 
As this is not realizable due to the growing size of the data packet, we seek to find a recursive 
estimation algorithm that gives the same accuracy as E Xf.(k)\, i.e., an algorithm that provides 

the same results as a central Kalman filter. Note, that in the given scenario, E|^ fc X^(A;)| is 

identical to E ja: fc Zfc(fc) j as the it's in (JS]) are known. 

In the following, an extension of the so called Hypothesizing Distributed Kalman Filter (HKF) 
is derived that calculates E |ai fc Zfc(/c) j based on locally pre-processed measurements when an 
assumption about the global measurement model has been met by the estimates. 

4-1- Hypothesizing Distributed Kalman Filter 

The idea of the HKF is to process local estimates in a transformed state space and filter 
measurements according to gains that are optimized according to a global measurement model. 
As the globalization of gains leads to biased local estimates in general, a correction matrix is 
maintained that allows to eliminate the induced bias. 

A detailed derivation of the algorithm including consistent mean-squared-error (MSE) matrix 
bounds is give n in [12]]. More precisely, the correctness of a central version of the HKF has been 
derived in [lol] and has been extended in [T2] to the concept of maintaining local variables that 
allow the calculation of a combined correction matrix. In this paper we limit ourselves to the 
presentation of key formulas and derive the subsequent incorporation of control inputs in more 
detail afterwards. All derivations and conclusions are given for the time-invariant system from 
Sec. [21 but are applicable one-to-one to the time-variant case. 

Let Gf, C {1, . . . , M} contain indices of sensors whose estimates are available at the controller. 
In order to apply the HKF, a hypothesis about the global measurement model (HGMM), i.e., 

(pi)- 1 « £ (cYi&y'a (ii) 

is necessary, which is chosen according to experimental data or is updated iteratively. Although 
an unbiased estimate is provided by the HKF in any case, the estimation quality in terms of the 
MSE matrix depends on the choice of the HGMM and is only equivalent to the globally optimal 
one when the HGMM meets the sum of the measurement models in (|11|) . 

We consider sensors that have no prior information and that initialize local values x l k and P^ 
by help of the first locally obtained measurement according to 

= L\y^ and P 1 = Pf with L| = P k (C i ) T (Qr 1 )' 1 . (12) 
When initial estimates ~Xq, P^ are given on sensor side, we set 

4 = Po^r 1 !^, and P = 

The variables are predicted according to 

4+i = and P fc+1 = AP fe (A) T + S* 
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Using L| from (|12|) and = Pfc(Pfc|fc-i) > the filter operation is given by 

xi=K fc 4 |fe _ 1+ 4|4 and P fc = ((P fc | fc _ 1 )" 1 +(PD" 1 )" 1 - 
For further processing, it is sufficient to transmit to the controller, where the fused estimate 



' - ^ Si (13) 



is obtained. It has been shown that xl equals the (central) linear minimum MSE result when 
equation (|lip holds exactly with G{ representing the sources of the estimates that are available [l(| • 
Otherwise, x{ is biased and must be corrected. 

To this end, a further variable A^ 1 is maintained that allows to reconstruct unbiased estimates 
independently of Gl and of the HGMM. The idea is to make sure that E{x k } = E{(A^)- 1 x|} holds 
for all local estimates at all time steps. Thus, we initialize the correction matrix with A^ 1 = ~L\C l 
when the estimates are initialized from measurements, or - when initial estimates are utilized - 
with A^* = Pi(P 1 ) _1 . The prediction of the correction matrix is led back to the last time step by 

-i 



and in the filter step, we set 



The correction matrix that corresponds to the fused estimate xj k is given by 

It is worth mentioning that A^* is not regular at the initialization step when the rank of C l 
is below the state dimension. However, this problem also occurs in the central processing schema 
and is likewise solved after fusing values that together cover the complete state. Apart from that, 
measurement models do not need to be known at remote sensors as potential differences between 
the HGMM and the actually utilized models in pip can be corrected by helps of A^. 

4-2. Subsequent Inclusion of Control Inputs 

Up to now, we have not considered control inputs from the state space model ( [!]) , which com- 
plicates the estimation process as the local sensors have to recursively estimate the state without 
having information about the deterministic inputs and thus, without knowing the complete system 
model. In the following, we determine the part of the control inputs that was not comprised in the 
processed measurements and needs to be added to the estimate. In a second step, we prove that 
the derived procedure provides globally optimal estimates iff the basic HKF is optimal. 

In order to simplify the calculations, we w.l.o.g. expect prediction and filter steps to be alter- 
nating. This assumption allows us to find sum formulas for the key variables. We define 

k— 1 T k 

G,.. fc = ( I] (A) T (K m ) T ) and (A,..^ 1 = f[ (A)' 1 , 
^ t=l ' t=l 
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with (A t ..j) 1 = I when t = lovt = l + l, and (A t ..j) 1 := A/ + i.. t _i when t > I + 1. Hence, the 
local estimation values and correction matrices are obtained as 

k k 

xi=J2G t .. k l4yi and A? =^G t .. fc LjC i (A t .. fc _ 1 ) _1 - 
The values of the fused estimate result in 

^ = E G ^(E%i) ( 15 ) 



and 



A*' =£G,.. fc ( 2 LjcfyiW.!)- 1 , (16) 



which verifies the key result from [12| that the estimation result of the HKF does not depend on 
whether information are processed locally or centrally as long as they are finally combined. The 
true state is given by 

fc-i 

^-k = E ( A t+i--fc-i B Mt + m) + A Q ..jfe_ 1 « (17) 
t=o 

and, hence, measurements can be represented as 

,t-l V 

= C * ( E (Ai +1 ..t_iBa, + uzj) + Ao-t-iSo ) +«j . (18) 
^ i=o ' 

By means of these formulas, we derive the part of the control inputs that is already included in xi 
and thus, determine the vector that still needs to be added in order to yield an unbiased estimate. 
We obtain: 

Theorem 4-1 The variable x{ = (A^) (xi + xt*) with 

^ = Y,G t .. k ^ {A)- l Bu t (19) 
t=o 



is an unbiased estimate of a system with dynamics ([T]) . 



Proof. Appendix A □ 



With xj from Theorem 14.1} we are able to reconstruct an unbiased estimate at the controller 
even if measurements are processed locally at the sensors. But, as we do not expect to have the 
same estimates available at every time step, we need the following corollary that allows to maintain 
the control input parts for every node separately. 
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Corollary 4.1 The varia ble x£ from JHJ equals the sum of node-specific variables 



fe-1 



# = £G t .. fe A?'(A) 'Bat. (20) 



t=0 

Proof. Simple matrix algebra using (|16p . □ 



It is worth mentioning that (|20|) depends on the control input as well as on the actually utilized 
model at the node (implicitly by A^ 1 ) and thus, both types of information have to be available at 
either the controller or the sensors. In this paper, we focus on the first case, where A^ 1 is obtained 
at the controller. Nevertheless, it would also be possible to communicate the control inputs directly 
from the actuator(s) to the sensors and transmit in combination with x k to the controller. 

Before we finish the derivations, we provide a Lemma concerning the optimality of the proposed 
procedure. 

Lemma 4-1 The HKF with subsequent control input inclusion is globally optimal when the HGMM 
has met the sum of the actually utilized measurement models from (|lip at all time steps. 

Proof. When the global measurement model is met, A^ = I, Vt € 1, . . . , k holds, and therefore, 
the fused estimate x k is obtained by 



lK f T\4+4 f )=d+4 f = E ^i) + EG,.. fc (A)- 1 B« i 

k 



By splitting up Kj and LJ, we obtain 

k 



E G *"* P * ( E ( C T ~ X tt + (P^-i) _1 Bn-, 



Mt-i , 



* =1 iec{ 



which equals the recursive information form of the Kalman Filter when all measurements are pro- 
cessed centrally and therefore, is the linear optimal solution. □ 



4-3. Application 

In the following, the derived approach is applied to the plant that has been depicted in Sec. [2j 
In order to minimize the communication and computation at the sensors, we calculate only x l k at 
the distributed nodes and maintain aft* and A^ at the controller. As estimates can be delayed or 
lost, the controller stores the measurement models and control inputs since the last time step of 
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successfully received variables g±j., t < k. By helps of this information, the estimation variables are 
predicted according to 

4 = G t .. k xt , A| 4 ' = G t .. fc A^(A t .. fc _ 1 ) _1 

and 

fc-l fc-i 



i=t 

k-l 



l=t l=t 

k-l 



I=t 

when no estimate of the current time step is available. If it is likely that transmissions fail, it is 
meaningful to store these predicted variables as future predictions are based on them. Otherwise, 
it is sufficient to hold variables A^,x" l ,Pt for the latest estimate for each sensor, and u k for 
each time step for which there is a chance that an estimate is predicted from. 

Although several matrices have to be stored and calculated at the controller, this is unlikely 
to be a problem in real-world applications where it is necessary to spatially separate controller, 
actuators, and sensors at the costs of a poor communication. Apart from that, optimal out-of- 
sequence algorithms, which are the alternative to the proposed algorithm, have at least comparable 
costs. 

It remains to discuss how the HGMM is chosen best. In order to apply the HKF it is possible to 
set the HGMM to an arbitrary matrix. However, as we are interested in optimizing the estimation 
quality at the controller, two strategies are meaningful in the context of the presented scenario. 
When the HGMM is set to the sum of the measurement models of all available nodes, i.e., = 
{1, . . . , M} in (jlip . the estimate at the controller is globally optimal when all estimates from one 
time step are received. This is especially meaningful when packet losses and delays are unlikely. 
Alternatively, the HGMM is chosen according to the expected "measurement quality", e.g., as the 
half of the sum of all measurement models when the sensor measurement models are similar and 
it is expected that slightly more than half of the transmissions fail. 

In plants with slowly changing communication quality, a more sophisticated approach can be 
employed. Although the exact procedure is out of the scope of this paper, it is reasonable to adapt 
the HGMM when the fused correction matrix (|14|) differs substantially from the identity matrix 
over multiple time steps. In this case, the HGMM should be reduced in dimensions with negative 
deviation to the identity matrix and increased otherwise. 

Independent of the specific strategy, the estimation result is almost optimal when the HGMM 



approximately meets the actually utilized models. In 11], it has been shown for multiple systems 
and measurement models in a sensor network consisting of ten nodes that the MSE of the HKF 
does not exceed the one of the linear optimal solution by more than six percent when the actually 
utilized global measurement model from (|11|) does not deviate by more than 40 percent from the 
HGMM. 

In summary, we have proven that by means of x^ from (|2U|) . it is possible to reconstruct an 
unbiased estimate from recursively gained information that is globally optimal when the HGMM 
fits to the measurement models of the available estimates, even if the control inputs are not known 
to the sensors but are subsequently included at the controller. In order to apply the proposed 
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algorithm, we have discussed strategies to maintain the relevant variables and to choose the HGMM 
appropriately. 



5. Conclusion 

In this contribution, a unified approach to sequence-based control for networked control systems 
was proposed. We have extended recent results on the optimal sequence-based control for systems 
with only one sensor sending raw data to the case of multiple sensors, that process measurement 
data in a sequence-based information framework. To this end, we introduced and extended the 
Hypothesizing Distributed Kalman Filter as distributed estimation algorithm. Although we have 
not presented an estimator that is equivalent to the central Kalman Filter in all cases, we have 
suggested methods for choosing a meaningful hypothesis about the global measurement model 
so that the estimator provides slightly suboptimal to optimal estimates when the communication 
model does not change abruptly. 

Future research will focus on relaxing the TCP-like assumption on the controller-actuator- 
channel and on suitable approximations of the multiple actuator scenario. The estimation procedure 
can be improved by handling measurement failures. Apart from that, we see great potential in 
applying the proposed algorithm to a real plant. 



Appendix A. Proof of Theorem 14.11 

Proof. In order to prove the theorem it is sufficient to show the equality between e|x{ + 
and e/a^'x^j. First, we note that the inner terms of x^J from (|16|) are given by 



G^Af'(A)- 1 i^G,,( £ LjC^A,,)- 1 . (A.1) 



The expected value of the estimate without considering control inputs is given by 

E{^} # EM E L * E fe}) f EG t .,ELiC i (^A w .,_ 1 B a + A„.,_ 1 Efe} 
which is transformed in order to factorize the control inputs to 

E ( E G t-k( E LjC^A^.^W + ^G^ £ LfGAo-.t-iEteo} . (A.2) 

1=0 v t=m ieG / J t=i igG / 

The inner term of (fAT2l) can be combined with the corresponding inner term of X/J that is given 
by G^Af^A)- 1 as 

E <M E LjG)^.,^ + G i+1 .. fc A^(A)- 1 <P £c t .. fe ( £ L'C^A,,)- 1 , 

t=t+i i&G { t=i ieG t 
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and thus, we obtain for e|x{ + j 

E (E G ^( E LjC^A^W + £ L;C*Ao.. t -iE{<So} . (A.3) 

With ([To]) and (fT7|) . the second term e|a^» a | is given by 

■ k \ , k—1 

Y, G t-k{ E LjC i )(A t .. fe _ 1 )~ 1 )(^A z+1 .. fc _ 1 B^ + Ao.. fc -iE{x } 



t=l ieG / ' M=o 



which is simplified with (A t ..k-i) 1 A/ + i..fc_i = (A^..j) 1 to (|A.3j) . □ 



Acknowledgments 

This work was partially supported by the German Research Foundation (DFG) within the 
Research Training Group GRK 1194 "Self-organizing Sensor- Actuator-Networks" and within the 
Priority Program 1305 "Control Theory of Digitally Networked Dynamical Systems". 

References 

[1] J. Hespanha, P. Naghshtabrizi, and Y. Xu, "A Survey of Recent Results in Networked Control Systems," 

Proceedings of the IEEE, vol. 95, no. 1, pp. 138-162, 2007. 
[2] A. Hekler, J. Fischer, and U. D. Hanebeck, "Sequence-Based Control for Networked Control Systems Based on 

Virtual Control Inputs," in To appear in Proceedings of the Control and Decision Conference 2012, 2012. 
[3] G. Liu, "Predictive Controller Design of Networked Systems with Communication Delays and Data Loss," IEEE 

Transactions on Circuits and Systems II: Express Briefs, vol. 57, no. 6, pp. 481-485, 2010. 
[4] L. Griine, J. Pannek, and K. Worthmann, "A Prediction Based Control Scheme for Networked Systems with 

Delays and Packet Dropouts," in Proceedings of the IEEE Conference on Decision and Control held jointly with 

the Chinese Control Conference, 2009, pp. 537-542. 
[5] J. Fischer, A. Hekler, M. Dolgov, and U. D. Hanebeck, "Optimal sequence-Based LQG Control over TCP-like 

Networks Subject to Random Transmission Delays and Packet Losses," Arxiv preprint, arXiv: 121 1.3020, 2012. 
[6] V. Gupta, B. Sinopoli, S. Adlakha, A. Goldsmith, and R. Murray, "Receding Horizon Networked Control," in 

Proc. Allerton Conf. Commun., Control Comput, 2006. 
[7] K. C. Chang, R. K. Saha, and Y. Bar-Shalom, "On Optimal Track-to- Track Fusion," IEEE Transactions on 

Aerospace and Electronic Systems, vol. 33, pp. 1271-1276, 1997. 
[8] W. Koch, "On Optimal Distributed Kalman Filtering and Retrodiction at Arbitrary Communication Rates for 

Maneuvering Targets," in Proceedings of the 2008 IEEE International Conference on Multisensor Fusion and 

Integration for Intelligent Systems (MFI 2008), 2008, pp. 457-462. 
[9] F. Govaers and W. Koch, "On the Globalized Likelihood Function for Exact Track-to- Track Fusion at Arbitrary 

Instants of Time," in Proceedings of the 14th International Conference on Information Fusion (FUSION 2011), 

2011, pp. 1-5. 

[10] M. Reinhardt, B. Noack, and U. D. Hanebeck, "On Optimal Distributed Kalman Filtering in Non-ideal Situa- 
tions," in Proceedings of the 15th International Conference on Information Fusion (FUSION 2012), 2012. 

[11] - — , "Decentralized Control Based on Globally Optimal Estimation," in Proceedings of the 51st IEEE Confer- 
ence on Decision and Control (CDC 2012), 2012. 

[12] - — , "The Hypothesizing Distributed Kalman Filter," in Proceedings of the 2012 IEEE International Conference 
on Multisensor Fusion and Integration for Intelligent Systems (MFI 2012), 2012. 



13 



